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>r^ ' We revisit a simple lattice model of aggregation in which masses diffuse and coalesce upon 
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contact with rate 1 and every nonzero mass chips off a single unit of mass to a randomly chosen 

neighbour with rate w. The dynamics conserves the average mass density p and in the stationary 
^1^' state the system undergoes a nonequilibrium phase transition in the (p—w) plane across a critical line 

W , Pc{u)). In this paper, we show analytically that in arbitrary spatial dimensions, Pc(w) ~ ^/w + 1 — 1 

exactly and hence, remarkably, independent of dimension. We also provide direct and indirect 
r~~« ' numerical evidence that strongly suggest that the mean field asymptotic answer for the single site 

mass distribution function and the associated critical exponents are super-universal, i.e., independent 
_^ of dimension. 
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^ ; I. INTRODUCTION 

^ . 

Nonequilibrium phase transitions [|l| occur in various systems including heterogeneous catalysis pL chemical reaction 

J5 ■ models [^, polynuclear growth models y, monomer-dimer models J 5|, models of fungal growth ||6[, non-equilibrium 

kinetic Ising models [0 , and branching annihilating random walks |S ] . A common feature in all the above is that the 

transition is from a state that has no activity to one that has continued activity. Such transitions are well characterized 

^ , by the critical exponents of Directed Percolation, Parity Conserving, or the DP2 universality classes. There are other 

Q ' nonequilibrium models that undergo phase transitions which do not belong to the above universality classes. These 

O , include boundary driven phase transitions [Bl and models whose steady states undergo a 'de-pinning' or 'unbinding' 
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transitions |10-1^. Recently a simple lattice model where masses diffuse, aggregate on contact, and also chips off a 



single unit of mass was studied |14 15f|. This 'chipping' model (CM) exhibits a nonequilibrium phase transition mA 15 



^ " in its steady state from a phase in which an infinite aggregate is present to one that has none. This nonequilibrium 

^— ^ , phase transition is in a completely different universality class compared to the other models studied in the literature 

, ' and mentioned above. The mathematical mechanism giving rise to the formation of the infinite aggregate at the onset 

^v , of phase transition was found to be very similar to that of the equilibrium Bose-Einstein condensation in an ideal 

f^ Bose gas. The difference is that in CM, the infinite aggregate or the condensate forms in real space as opposed to the 

^D ! Bose gas where the condensation takes place in the momentum space. Besides, the phase transition in the CM occurs 

^^ ' even in 1-dimension as opposed to the Bose gas where the condensation occurs in 2 and higher dimensions. A slightly 

different off-lattice version of the CM was studied earlier within rate equation approach in the context of aggregation 

in dry environments |lq] . A directed version of the CM where masses move asymmetrically only along one direction 

was studied in Ref. |lq| and its critical properties were found to belong to a different university class from that of the 

'^ undirected CM. This directed CM also appeared recently in the context of a traffic model with passing [[l7[ . 

^ The undirected CM is defined on a d-dimensional hyper cubic lattice with periodic boundary conditions. Starting 

O ' from a random distribution of nonnegative integer masses at each lattice site, the system evolves in continuous time 

via the following microscopic processes. In an infinitesimal time interval Ai, (i) with probability At, the mass at 

each site hops to one of its neighboring sites, chosen at random, (ii) with probability wAt, a unit mass is chipped off 

from an already existing mass at each site and added to one of the neighboring sites, again chosen at random, and 

5J] ■ (iii) with probability (1 — (1 + w)At), each mass stays at its original site. Following the steps (i)-(iii), the masses at 

5^ \ any given site add up. Since this system is closed to the environment (periodic boundary condition), the total mass 

is conserved by the dynamics. Thus, there are two parameters in the problem, the average mass per site p, and the 

ratio of the chipping rate to the rate of hopping as a whole, w. In the long time limit, the system evolves into a time 

independent steady state. The steady state single site mass distribution function P{m), i.e., the probability that a 

site has mass m when i — *■ cx), was shown to undergo a phase transition in the p-w plane |l4,Qq]. There is a critical 

line Pc(w) in the p-w plane that separates two types of asymptotic behaviors of P{m). For fixed w, as p is varied 

across the critical value Pc{w), the large m behaviour of P{m) was found to be, 



3-"/™* p < pc{w), 

P{m)^{m-- p = PcH, (1) 

^^'^ + infinite aggregate p > Pc{w). 

Thus, the tail of the mass distribution changes from having an exponential decay to an algebraic decay as p approaches 
Pc from below. As one increases p beyond Pc, this asymptotic algebraic part of the critical distribution remains 
unchanged but in addition an infinite aggregate forms. This means that all the additional mass (p — Pc)y (where V 
is the volume of the system) condenses onto a single site and docs not disturb the background critical distribution. 
This is analogous, in spirit, to the condensation of a macroscopic number of bosons onto the single k = Q mode in an 
ideal Bose gas as the temperature goes below a certain critical value. This infinite tower, i.e., the single site carrying 
macroscopically large (proportional to volume) mass is mobile. Mathematically this means that for p > Pc{w), 
P{m) ^ m~'^ + y (5(771 — aV) where a = p — pc- This last delta-function peak contributes only to order 1/V in the 
integral J P{m)dm indicating that the aggregate is contained in thermodynamically few sites (indeed a single site) 
but contributes a finite amount a in the integral J mP{m)dm even in the thermodynamic limit V — > oo. 

These results were found both analytically within a mean field approximation which involved ignoring all correlations 
between masses, and numerically in one dimension. Within the mean field approximation, the locus of the critical 
line was found ||l4|] to be, Pc{w) = VuT+T — 1 and the exponent r = 5/2 |14|Jl6|] . In one dimension, the numerically 
obtained critical line p4| was found to be close to the mean field phase curve. Even the exponent r, determined from 
a simple linear fit on the log-log plot of data from relatively smaller size lattices, was found to be t « 2.33, rather 
close to the mean field exponent 2.5. This raises the question whether or not the mean field answers for asymptotic 
behaviors of P{m) are exact even in one dimension. On the other hand, one can show explicitly (see section II) that 
there exist nonzero correlations between masses in any finite dimension, even in the thermodynamic limit. Thus, one 
is confronted with a puzzle. 

The purpose of this paper is to shed new light on this puzzling issue. We first prove analytically a remarkable result 
that the mean field phase boundary, Pc{w) = \/w + 1 — 1, is indeed exact and independent of the spatial dimension 
d. This, of course, still does not prove, but hints, that the exponent r may also be independent of d. However, we 
provide rather unambiguous numerical evidence, in conjunction with several direct and indirect checks, which suggest 
that even the exponent r = 5/2 is super-universal and independent of d. Thus, our results seem to suggest strongly 
that even though there are nonzero correlations between masses in finite dimensions, these correlations do not affect 
the asymptotic behavior of the single site steady state mass distribution P{m). However, the possibility remains that 
other higher order correlation functions, such as the joint distribution of two masses P{mi,m2) will depend on the 
spatial dimension d. 

The paper is organized as follows. In section II, we show analytically that in arbitrary dimensions d, the locus of the 
phase boundary is independent of d and is given by the mean field expression, Pc{w) = -y/uT+T— 1. In section-Ill, we 
re derive analytically the mean field expression for the mass distribution P{m) by a method different from that used 
in |]l4| and compare these analytical expression with the numerical results obtained in one and two dimensions. In 
section IV, we do a finite size scaling analysis that provides additional evidence that the exponent r is super-universal. 
Finally, we conclude with a summary and discussion in section V. 

II. EXACT PHASE DIAGRAM IN ARBITRARY DIMENSIONS 

As mentioned in the introduction, the steady state of the CM undergoes a phase transition in its parameter space 
{p-w plane) across a critical line Pc{w) in all dimensions. In this section, we compute Pc{w) exactly in arbitrary 
dimensions by analyzing the two point equal time mass-mass correlation function, C(x, t) = (to(x, t)?7i(0,i)). Let 
77(x, x', t) denote the mass transferred from a site x to a neighboring site x' in the time interval between t and t + At. 
Clearly, r/(x, x',t) is a random variable that takes the following values, 

{to(x), with prob. Jd^^'' 

1 - '5m(x),o, with prob. j^wAt, (2) 

0, with prob. (l-i±f Ai), 

where 2d is the number of neighbors of any given site. The Kronecker delta function in the second line of Eq. (g) 
indicates that a chipping of a unit mass can take place provided the mass m(x) is a positive integer bigger than 0. 
Then, the evolution of mass at site x can be written as, 

m(x,t + At) = m(x,t)-^77(x,x',t)-f ^77(x',x,t), (3) 



where the sum is over the neighbors x' of the site x. The second term on the right hand side of Eq. (g) describes the 
outflow of mass from x while the third term accounts for the inflow of mass into x from neighboring sites. It is quite 
straightforward to write down the 2 point correlator for rj to order At. Suppressing the explicit t dependence in 77 we 
find, 

(77(xi,x'i)?7(x2,x'2)) = — (m^(xi) + w{l - 5,„(xi),o)) (^xi.xs^x'i.x's At. (4) 

Using Eqs. (0) and (0), the evolution equations for the two-point equal time correlation function, C{x,t) = 
(77T,(x, t)m(0, t)) can be written down. Multiplying Eq. (0) by m(0,t + dt) and taking average on both sides, and 
putting all time derivatives to zero in the steady state, we get, 



x' \ x' / 



-(C(0)+«;s)Ux,o-^}_^<5x,xol, (5) 

where Xq denotes the neighbors of the site 0. Also, D(x, i) = (^7i(x, i)(5„i(o_t)^o)j and s = 1 — ((5m(o),o)j is the 
probability of a site having non zero mass. Thus, in the CM, the two point correlation functions do not form a closed 
set of equations making it difficult to solve for C(x, t) exactly. This is unlike many other models where two point 
correlations do form a closed set of equations and hence are solvable. A few examples include the 1-d Glauber model 
|l8| , asymmetric random average process ||l^,g2li the Takayasu model of aggregation [gl| and the q model of force 
fluctuations in bead packs pi|,p2| . 

Remarkably, however, Eq. (g|) allows for the solution, 

C(x) == w[D{x.) - s] , X 7^ 0. (6) 

This solution is also the unique solution. To see this, observe that the homogeneous part of Eq. M) is the Laplace's 
equation V^(C(x) — wD{x)) — 0, with the boundary condition that (C(x) — wD{x.)) is a constant as | x |— > 00. Since 
the solution Eq. (ra) satisfles the inhomogeneous part too, as well as the boundary conditions, it is the unique solution. 
Note that the above solution Eq. (||) is also valid on a finite lattice. Summing Eq. (g) over all x 7^ and using the 
fact that the conserved total mass is given by, X]x"^(^) ~ P^ (where V = L'^ is the volume of the system), we get 
the following exact equation, 

P' - ^ - wp{l -s)- ws. (7) 

This equation is reminiscent of Bose-Einstein condensation in ideal Bose gas. In the low density phase, we expect that 
the system reaches a stationary state in which (tti^) is a finite number of order 0(1). Therefore, in the thermodynamic 
limit 1/ ^ 00, the second term on the left hand side of Eq. (0) drops out and we get, 

p^ — wp{l ~ s) — ws. (8) 

Note that the Eq. (0) could have been obtained from Eq. (pi) if we had assumed that the two point correlation 
functions decouple, i.e., C(x) — {■m{0)m{x.)) — (TO(0))(m(x)) — p^ and -D(x) = (r7i(x)(5„j(o),o) = /'(I — s)- However, 
there is apriori no reason for the two point correlations to decouple. Our derivation of Eq. (|^) does not rely on this 
decoupling. 

From Eq. (g) we get, 

w(p+ 1) 

According to Eq. (O), as one increases the density p keeping w fixed, the occupation probability s first increases with 
p, attains a maximum at p = \/w + 1 — 1 (obtained by setting 4^ = in Eq. (O)) and then starts decreasing with 
increasing p. However, it is clear that s, the probability that a site has nonzero mass, must be a monotonically non 
decreasing function of p. Hence we conclude that Eq. (m is valid as long as p < pc = ^/w + 1 — 1. For p > pc, the 
basic assumption (rr?) jV ^- as F ^ 00 breaks down and Eq. (pf) ceases to be valid. 
Thus, the critical density is given by, 



PeH = \/rT^-i, (10) 



and remarkably, it is independent of d and not surprisingly, therefore, coincides with the mean field expression |14| , |15[ . 
For p < Pc, s is given by Eq. (pl). As p increases from to Pc{w) (for fixed w), s increases monotonically according 
to Eq. (P) up to the value Sc given by. 



For p > Pc{w), s does not increase any further and sticks to its value Sc- Putting s = Sc in Eq. m) and using the 
expression of Sc from Eq. (pl|), we get for p > pdw), 

lim(^) = (p-p,)2. (12) 

Thus, for p > Pc, {m?) becomes macroscopic, i.e, proportional to volume. Since s, the fraction of occupied sites, does 
not increase anymore for p > pc, this indicates that all the extra mass (p — Pc)V condenses onto a thermodynamically 
negligible number of sites (indeed, a single site only) with density ^ ■^, leading to the macroscopic behavior of 
(m^) ~ y[ip~Pc)V]'^ ^ [p — pcYV . This is similar in spirit, though not in details, to the Bose-Einstein condensation 
where below a certain temperature, the number of particles in the fc = mode also become macroscopic. 

Let us conclude this section by stressing on an important point. We note that the mean field solution (assuming 
decoupling) for the stationary two point correlation function, C(x) = p'^ for x 7^ and -D(x) — p{\ — s) for x 7^ with 
s satisfying Eq. (0), is indeed an exact solution of Eqs. (||) and (||). However, this need not be the only stationary 
solution. Besides, even if the mean field answer for the two point stationary correlation function is the correct one, 
it still does not prove that the mean field theory is exact. For example, one can show [g^ that indeed 3 and higher 
point correlation functions do not decouple. In any case, the main result of this section, namely the derivation of the 
exact phase boundary, does not rely on whether the correlation functions decouple or not. 

III. COMPARISON WITH MEAN FIELD THEORY 

In the previous section, we proved that the mean field phase diagram is exact in any dimension. This, of course, 
does not prove but suggests that, perhaps, even the mean field expression for the distribution P{m) may also be 
asymptotically exact in all dimensions. In this section, we try to provide evidence in favor of this hypothesis. For this 
we first re derive the mean field expression for P{m) and compare it with the numerical results obtained in 1 and 2 
dimensions. 

In Ref. |14| , the steady state single site mass distribution function P{m) was computed analytically by assuming 
that the joint distribution P{mi, 7712), the probability that two consecutive sites have masses mi and 7712 respectively, 
factorises, i.e., P(rni, 7772) = P(7n.i)P(7n.2). With this assumption, P{m) was shown to satisfy a closed set of equations 
which were then solved via generating function method yielding results mentioned in Eq. (n^) with t = 5/2 and 
Pc{w) = y/w + 1 — 1. In this section, we first re derive the mean field results by a different method that requires lesser 
restrictions than the product measure used in Ref. ||l4| ]. 

Here we use a technique used before for solving the mass distribution function in other models of aggregation |2^] 
as well as the w — Q limit of the CM |2J]. We consider the CM on a 1-d lattice. Let P{mi, 7772, . . . , 777„) denote the 
joint probability that 77, consecutive sites on the lattice have masses mi, 7712, • ■ •, m,n respectively in the stationary 
state. We define two generating functions. 
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^«= E ••• E a;™^+-+""^(™i,---,"^n), (13) 

00 00 

^-= E ••■ E a;™^+-+""^(™i,---,™n,0). (14) 
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Here Z„ — ^x™i +■+'"") is an unconditional average but F„ = ^^rni+...+rn„'^^ jg ^ conditional average where the 
(tt- + l)-th site is conditioned to have mass. Using the dynamics of Trie's (as given by Eq. (H)) and following steps 
similar to those used in Refs. |23J24], one can write down the evolution equations for ZnS. In the steady state, when 
all time derivatives go to zero, we get, after some algebra. 



Z„+i - 2Z„ + Z„_i + w ( ^^^Z„ - ^-^Yn-i + (1 - a;)i;0 ) = 0, (15) 



with the boundary conditions, Zq — 1 and Yq = -P(O) = 1 — s, the probabihty of having no mass at any given site. 

If we now make the assumption that P{rni + . . . + to„ — m^Q) = P{mi + . . . + m„ ~ m)P(O), i.e., Y^ = P{0)Zn, 
then we get equations that contain only the Z„^s. As mentioned before, this assumption is less strict than the product 
measure approximation as it requires the factorization of only a special conditional probability where a site at the 
beginning of a string is empty. In order to determine P{m), we need to compute Zi{x) which, by definition in Eq. 
(Q) is the generating function for P(m)'s, 

oo 

Zi{x) = Yl x"'Pim)- (16) 

Since Zi depends on other Z„'s, we need to solve the full Eq. (^5|) for all n. Eq. (p^ ) can be solved by standard 
generating function method. Let G{x,y) — X^^i ^n{x)y"- Multiplying Eq. ( |l5| ) by y" and summing over n, we get, 
after straightforward algebra, 

(^/ X ^ y [wP{0)yil -x)~xy + xZ^{x)\ 

x{\ — y)^ + wy{l — x)^ + wP{0)y{l — x){x — y) 

For a fixed x, when considered as a function of y only, G{x, y) has two poles. 



— 2x + w{l — x){l — Sx) ± (1 — x)y/wy/w{l — SX)'^ — isX 

y'''= 2ai~s)wii~x)-x) ■ ^^^> 

For a fixed w and s, | y2 |< 1 for small values of x. This implies that Z„ ^| 2/2 | " for large n. However we cannot 
have a diverging probability for large n. Hence this pole must be canceled off by the numerator of G. This 'pole 
canceling' mechanism was also useful in deriving exact results in other recently studied aggregation models po| , pi| . 
Demanding G{x, 2/2) = 0, we get the following expression for Zi(x), 



2x — w(l — a:;)(l — sx) + (1 — x)y/w\/w{l — sxY — 4sx 
Zi[x) — - . (19) 

The coefficient of x"* on the right hand side of Eq. ( pj| ) will then give the desired distribution P{m). 

The expression for the generating function Zi{x) = X]m=o ■^™^("^) ^^^ ^^- ( M ) is identical to the one derived in 
|l4| using the approximation of product measure. Thus, the two methods, though different in details, yield the same 
P{m) for any m and not just for large m. Thus, this result for mean field P{m) seems to be extremely robust and 
does not depend on the details of how the mean field assumption is incorporated. 

The asymptotic properties of P[m) for large m can be derived by analyzing the behaviour of Zi{x) near a: = 1 |lj] 
and one recovers t_he results in Eq. (|l|) with r = 5/2 and Pc{w) = \/w + 1 — 1. However, by expanding the expression 



for Zi{x) in Eq. (g^ in powers of x using Mathematica, we have determined P{m) for all m. In the aggregate phase 
(p > pc), we set s = Sc in Eq. ( p^ ) and calculate the distribution P{m) by expanding in powers of x. In Fig. (1), 
we compare this analytical mean field answer for all m with the numerical results obtained in 1 and 2 dimensions in 
the aggregate phase (p > pc)- Note that Fig. (1) shows only the power law part of the spectrum. The numerical 
data is for V = 900 in 1-dimension and for a 30X30 lattice in 2-dimensions. The two curves for 1 and 2-d are almost 
indistinguishable from each other. While the numerical data matches excellently with the mean field answer for small 
771, it has a small deviation for larger masses. This deviation at large m is due to finite size effects. To confirm this, 
we also did simulations for larger lattice sizes up to V^ = 2000 in 1-dimension. For a fixed mass, we confirmed that 
the deviation from the mean field decreases with increasing V (see the panel inside Fig. (1)). We also compared the 
mass distribution in the exponential phase (p < pc) with the mean field prediction. Again excellent agreement is seen 
(see Fig. (2)). These results thus provide strong evidence that the mean field expression for P{rn) is exact in all 
dimensions and therefore the exponent t = 5/2 is also super-universal, i.e., independent of d. 

IV. FINITE SIZE SCALING AND INDIRECT NUMERICAL CHECKS 

In this section, we provide further indirect numerical checks which again strongly suggests that the mean field 
answer for P{m) is indeed exact. We start by making a reasonable finite size scaling ansatz for P(m, V) (where V is 
the volume of the system) in the aggregate phase p > Pc for large 777., 



1 / 771 \ 1 

Pi-^^ y) - -ff (^) + ySim~{p^ p.)V) , (20) 

where the exponent is a crossover exponent and the delta function peaked at (p — Pc)V indicates the aggregate 
containing a macroscopic amount of mass. The power-law part of the mass distribution gets cut off at ~ V^ for 
finite V. Since the largest mass inside the power law part is much less compared to that in the aggregate, we get the 
inequality, < 1. 

Also, since the mass density, J mP{m)dm = p is finite. Substituting Eq. (EO) in this integral, it is evident that 
for large V, the integral will converge provided r > 2. Also, from Eq. (po|), one finds that the second moment, 

■i^^ « (p — pcY + 0{V'^''^~'^'^~'^). The first term is from the aggregate and the second term is from the power-law 
part. This is consistent with the exact result Eq. (n2[) provided, 0(3 — r) < 1. This provides an upper bound for r, 
i.e., T < 3. Thus, we have the exact bounds for r, 2 < r < 3. The mean field results as well as numerical simulations 
in 1 and 2-d are consistent with these bounds. 

Next we derive a scaling relation between r and (p. This is obtained by demanding that P(m, V) is normalized, 
J P{m,V)dm — 1. Substituting Eq. ( pO[ ) in this integral, we get, 

±fQ)dm = c-^, (21) 

where c is a constant of 0(1), and mi ^ 0(1) , is the mass beyond which the scaling starts holding. Differentiating 
Eq. ( |2l| ) twice with respect to V and making a change of variable, we arrive at the relation, 

f"(y)y'-^dv = -{l-<f>)V*--'~\ (22) 

where f"{y) = -;jt-- The reason for twice differentiation is as follows. Since 2 < r < 3, one can replace the lower limit 
of the integral in Eq. (E2h by in the V ^ oo limit and there is no divergence from the lower cut off. Once the lower 
cut off is replaced by 0, the integral on the left hand side of Eq. ( p2| ) is of 0(1). Comparing this with the right hand 
side of Eq. (|2^), we immediately arrive at the scaling relation, 

0(r-l) = l. (23) 

If T = 5/2, Eq. (H) would indicate = 2/3. 

We found that the cleanest way to measure the exponent <j) is via the following indirect finite size method. On a 
finite lattice, the critical value of the occupation probability of a site, sdV) will differ from the y = cx) value, Sc (oo). 
One can assume a reasonable finite size correction of the form, 

'^^y)-'T~^, (24) 

where a is a constant and is a new exponent. Using this ansatz along with the expression of (m^) obtained from 
Eq. ( pol ) in Eq. (0), we find y^(3-T)-i ^ y-^ ^ implying another scaling relation, 

6( = l-0(3-r). (25) 

If T =: 5/2 and cj) = 2/3, we get from Eq. (p5|), 9 — 2/3. In Fig. (3), we plot, in one dimension, Sc{V) for various 
values of V and indeed find that 6 = 2/3. This provides an indirect numerical check on these scaling assumptions as 
well as suggests, once again, that r = 5/2 is a super-universal exponent. 

Another indirect check can be done by mapping the one dimensional CM to a fluctuating interface model |1^ . We 
outline the procedure in brief. The flrst step is to map the CM on to a system of hard core particles moving on a 
ring. For this, interpret rrii, the mass at site i, to be the gap between the (j — 1)*'' and i*'' particle. Then the chipping 
move corresponds to a particle jumping forward one step provided the target site is empty, while the aggregation 
move corresponds to a particle making a long range jump to the site adjacent to the particle nearest to it. There 
is a standard procedure to map a lattice gas configuration to an interface configuration ||25| . Let Uj = 1(— 1) if a 
particle is present (absent) at site i. Then hi = ^j^inj. While the chipping move corresponds to local moves of 
the interface, the aggregation move corresponds to nonlocal moves of the interface. The width of this interface was 
monitored numerically p5[ as a function of time t. On the critical line Pc{w) of the CM, the width was found to have 
the scaling form p5|] , 



wiV,t)^V^f(^^j, (26) 

with X « 0.67 and z « 2.0. 

We show here that the roughness exponent x can be related to the exponents r and cj) of the CM. We first map the 
CM in 1-d to an interface model with the height field hi — J2]=i '^j- Then (hi) — ip and {hf) w i{m?) + i{i — l)p^. 
If we approximate {mirrij) by p^ (assuming that the mean field P{m) is exact), the width is then simply uP = 
1/F ^,-^j^ (ft,i — {hi)Y ~ V^i+"^(3-r) _ ]/2x^ This implies that provided {minij) « p^, one gets the scaling relation, 

X = ^[l + 0(3-r)]. (27) 

If T = 5/2 and ~ 2/3, we would get from Eq. (P7|), x = 2/3, in excellent agreement with the numerical value 
X ~ 0.67 |l^. This is further evidence for r = 5/2 and <f) = 2/3. 

V. SUMMARY AND OUTLOOK 

In this paper, we have studied a simple stochastic lattice model where masses diffuse as a whole and coalesce upon 
contact with rate 1 and every nonzero mass chips off a single unit of mass to a neighbour with rate w. The mass 
density p is conserved by the dynamics. This model undergoes a nonequilibrium phase transition in the (p — w) plane 
across a critical line Pc(w). We have shown analytically that pdw) = y/w + l — 1 in all dimensions. We also provided 
both direct and indirect numerical evidence that strongly suggest that the mean field answer for the single site mass 
distribution function P[m) might be exact in all dimensions and that the exponent r — 5/2 is super-universal. 

However we would like to stress one important point. Even though the single site distribution P{m) may be 
given exactly by the mean field answer, that does not prove that mean field theory or product measure is the exact 
stationary state in all dimensions. In this sense, the CM is different from the recently studied g-model of force 
fluctuations pq | where the product measure is exact in the stationary state. More precisely, in the g-model, the 
evolution equation for single site mass distribution P(m) involves the 2-point distribution P{jni,m2)- Similarly the 
equation for 2-point involves the 3-point function P(mi,m2,m3) and so on. However, if one assumes the product 
measure, i.e., P{mi,m2, ■ ■ -mn) — P{mi)P{m2) ■ ■ .P(?7i„) for all n, then it was shown ||26| that this ansatz satisfies 
all the equations of the hierarchy. In the CM, in a similar way, one can write down the full hierarchy of equations 
satisfied by the n-point distribution functions. However, unlike the (/-model, the product measure ansatz does not 
satisfy all the equations of this hierarchy ^^. Nevertheless, the expression of P{m) obtained from the first equation 
of this hierarchy (the equation that involves only P{m) and P(mi, TO2)) seems to be extremely close to the numerical 
answer. This suggests that in the stationary state of CM, correlations between masses appear only in the 3 or higher 
order correlations but seem to be absent at the 2-point level. This remarkable fact was also noticed recently in another 
aggregation model namely the asymmetric random average process with sequential updates |19| , |20|| , suggesting that 
such unusual correlations may be more generic and less exceptional. 

In this paper we have studied the undirected CM. As mentioned in the introduction, the directed CM also has a 
qualitatively similar phase transition in the steady state though the associated critical exponents are entirely different 
from the undirected one. Also, the phase boundary in the p—w plane of the directed CM was found be quite different 
from the mean field phase boundary urn. A proper understanding of the directed model remains an outstanding 
challenging problem. 

We thank M. Barma, D. Dhar and S. Krishnamurthy for useful discussions. 



10> 



^ 

K 



10-' 



10-2 



10-3 



10-* 



10-" 



10- 



A V = 400 

. V = 800 

n V=1000 

X V = 2000 

— mean field 




Li mean field 

V= 30x30; 2-dim 

V=900; 1-dim 



10 



100 



rri 



FIG. 1. The power law part of the steady state mass distributions P{m) in 1- and 2-diniensions are compared with the mean 
field answer in the aggregate phase p > pc- The data is for ui = 3 and p = 10. The critical value is pc = 1 for w = 3. In the 
inset panel, we show the convergence of the probability distribution to its mean field value as the system size is increased in 
one dimension. 
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